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^—N ' We study the role of small-size instantons in the determination of the topologi- 

Q«^ ■ cal susceptibility of the 2-d 0(3) a model on the lattice. In particular, we analyze 

C^ , how they affect the non-perturbative determination, by Monte Carlo techniques, of 

pi 1' the renormalizations on the lattice. As a result, we obtain a high-precision non- 

r^ • perturbative determination of the mixing with the unity operator, finding good 



X 



agreement with perturbative computations. We also obtain the size distribution 



H ' of instantons in the physical vacuum up to very small values of the size in physical 

units, without observing any ultraviolet cut-off. Moreover, we show by analytical 
calculation that the mixing of the topological susceptibility with the action density 
is a negligible part of the whole non-perturbative signal. 



I. INTRODUCTION 

Simulating a theory on the lattice is the only reliable instrument to study its non- 
perturbative aspects: the lattice theory, merely being an UV regularized version of the 
theory on the continuum, takes into account all non-perturbative fluctuations. Some 
quantities of physical interest, which can be extracted from the lattice, are related to 
vacuum expectation values of composite operators. Lattice Monte Carlo simulations 
give a numerical estimate of the cut-off-dependent bare expectation values, while the 
physical quantities on the continuum are the renormalized expectation values, which are 
cut-off-independent. So, the physical quantities can be determined from the lattice only 
if the renormalizations of the lattice-regularized operators are completely under control. 

Perturbation theory has been, up to now, the only mean for the evaluation of the 
lattice renormalization constants. In particular, perturbative techniques, combined with 
Monte Carlo simulations, have been applied to the problem of the determination from 
the lattice of the topological susceptibility of the QCD vacuum [|I|] . The situation 
has recently changed, since a new method for the determination of the renormalizations 
of the lattice topological susceptibility has been found [Q. This method, known in 
literature as the "heating" method, is fully non-perturbative since it relies only on Monte 
Carlo techniques. Applications to CP^^^ models [Q and SU{2) Yang-Mills theory [Q 
have shown a clear agreement with perturbative calculations 0. In SU(2) it has 
been possible to obtain also an indirect determination of the gluon condensate, finding 
agreement with previous standard Monte Carlo determinations [^. 

The 0(3) a model is the simplest model displaying a non-trivial topology, and so it 
appears as the best laboratory to test the heating method. The drawback is that this 



topology is suspected to be pathological: the semiclassical approximation [|IOl shows 
a small-size divergence in the instantons contribution to the partition function. Such 
an ultraviolet dominated topology is expected to strongly affect the heating method, 
which relies on a decoupling between short-range perturbative modes and topological 
fluctuations, assumed to be long ranged. 



In fact, the method had its first apphcation just in 0(3) a model |Q p, but the 
statistical fluctuations in the numerical results of those first works are likely to mask an 
eventual exotic behavior of the model; so, a more accurate investigation on the matter 
is asked for. 

In the present work we perform a careful analysis of the heating method in the 0(3) a 
model to check if, and to what extent, the ultraviolet feature of the topology spoils the 
non-perturbative determinations of the renormalizations of the lattice susceptibility. 

II. THE MODEL 

The 2-d 0(3) o model or OP^ model is described by the Lagrangian: 

>C = - d^4){x) ■ d^(f){x) , (1) 

where 0(x) is a three component real field satisfying the constraint (j) ■ (j) = 1. 

The 2-d 0(3) a model plays an important role in quantum field theory because it 
resembles in several aspects the 4-d non-Abelian gauge theories: asymptotic freedom, 
non-perturbative behavior in the infrared region with spontaneous mass generation, 
non-trivial topological structure. 

The topological charge of a spin field 0(x), Q, is the number of times 0(x) winds the 
sphere 5^. It can be expressed as the integral over the space-time of a local operator, 
Q(x): 



Q{x) is the divergence of a topological current K^ [0 [^ 



Q{x)=d,K,{x) . (3) 

All the classical solutions with non-trivial topology, the fc-instantons, have been 
explicitily found [l^ . At a quantum level, the only available prediction comes from the 



semiclassical approximation, which gives for the size distribution of instantons in the 
physical vacuum: 



Vdp ^ p3 



^-4.,^2M^ , (4) 



where M is the cutoff mass; renormahzation group theory implies 

/(x(pMr-*, (5) 

where /3o = 2^ is the one-loop coefficient of the P function; as a result: dN/dp oc 1/p. 
The generalization for CP^-i is: dN/dp oc p^^^ For SU{N), Eq. (§ and Eq. (g) 
become respectively 



dN ^ s.yg--sfipM) 

Vdp ^ p5 



^-8.V.^^-8M^ (6) 



and 



/ oc ipMY'-'^^ ; (7) 

it follows: dN/dp oc p^^-^. For both CP^-i, with A^ > 4, and SU{N), the size 
distribution is suppressed at small sizes with a power law, while in the 0(3) a it diverges. 
This makes 0(3) a model singular, and a behavior radically different from the theory of 
physical interest, QCD, is expected at small distances. 

III. THE FIELD THEORETICAL METHOD 

We regularize the theory on the lattice by taking the Symanzik tree-level improved 
action ill 



s^ = -/3y: 



X,fJ, 
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-0(x) ■ 0(x + p) - Y^0(a;) ■ 0(a; + 2p) 



In the field theoretical method, a lattice topological charge density operator is defined 
as a local operator having the appropriate classical continuum limit []T3|; our choice is 



Q^ = ——e^^eijk(t)i{x){(t)j{x + p) - (t)j{x- p)){(t)k{x + v) - ^^(x - z/)) . (9) 
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It has been shown in the framework of perturbation theory ||^ that the operator 
Q{x) defined on the continuum is invariant under the renormahzation group. 

A finite multiphcative renormahzation connects the matrix elements of Q^{x) with 
those of Q{x) defined on the continuum. The antisymmetry of Q^{x) forbids mixings 
with any other 0(3) invariant operator of dimension two or less. So, the connection 
with the continuum is: 

Q^(x) = a^Z{/3)Qix) + Oia") . (10) 

The topological susceptibility is defined on the continuum as the correlation at zero 
momentum of two topological charge density operators, Q{x) : 



X = J d'x{0\T[Q{x)Q{0)]\0) . (11) 



The prescription defining the product of operators in Eq. (|ri|) is [16 



(o| T[ g(x)g(o) ] |o) = d,{o\ t[ k^{x)q{q) ] |o) . (12) 

This prescription eliminates the contribution of possible contact terms (i.e. terms pro- 
portional to the 6 function or its derivatives) when x — *> 0. 
The lattice-regularized version of x is 

X^ = (EQ^(^)Q^(O)) = ;^((EQ^(^))') > (13) 

where L is the lattice size. 

A prescription equivalent to Eq. (|l^) does not exist on the lattice, and therefore the 
contribution of the contact terms must be isolated and subtracted. These contact terms 
appear as mixings with the action density S{x) and with the unity operator /, which 
are the only available operators with equal dimension or lower. In formulae 

X^(/?) = a^ Z{(5fx + a' A{(3) {S{x)) + P{[3) (/) + 0{a') , (14) 



where a is the lattice spacing. In Eq. (|Tj) the quantity {S{x)) is intended to be the 



non-perturbative part of the expectation value of the action density, i.e. it is a signal of 
dimension two. 



Z{P) and P{P) have been calculated perturbatively in Ref. up to the order l//3^ 
and l/P^, respectively. 

A relation similar to Eq. (|l^ holds in SU(2) with the appropriate powers of a and 
{S{x)) -^ {g'^ / {At^Y F'j^^F^^{x)) , the gluon condensate. In SU{2), the mixing with the 
gluon condensate is an appreciable portion of the whole non-perturbative signal in the 
scaling region, and can be detected through x^ H- 

We realize that the situation is radically different in the case of the 0(3) a model. 
Indeed, calculating the leading perturbative contribution to the mixing coefficient, and 
using a large N result for {S{x)) [|l7l (see Appendix for the details of the calculation 



and the discussion), we find that the mixing term is smaller than one thousandth of the 
total non-perturbative signal in the scaling region observed in Ref. [Q. We therefore 
argue that the contribution to x^ coming from the mixing term can be safely neglected. 

IV. NUMERICAL RESULTS 
A. Heating and cooling 

Now we want to give a brief account of the heating method (for Ref., see [Q [Q). 
This method allows a non-perturbative determination of the lattice renormalization 
constants, Z{(3) and P{(3), which have their origin in the short-range fiuctuations (/ ~ a). 

Ensembles of configurations {Ct} are constructed on the lattice, each configuration 
of the ensemble being obtained by performing a sequence of t local Monte Carlo sweeps 
starting from a discretized classical configuration Co - a large instanton or the fiat con- 
figuration. For small t, the heating process does nothing but "switch on" the small-range 
fiuctuations which are responsible for the renormalizations: when the starting configura- 
tion is a large instanton (fiat configuration), measuring Q^{x^) on the ensembles {Ct}, 
a plateau at the value of Z{f3) {P{f3)) is expected after a certain time, not depending 
on /?, corresponding to the time of thermalization of quantum fiuctuations. 

If the heating is protracted, the (local) algorithm of thermalization generates fiuctu- 
ations of ever increasing size according to the random walk law Z^ oc t. 



In the case of the heating of the flat configuration, when the ranges I ^ ^eq - the 
equihbrium correlation length - are reached, the contribution of the mixing with the 
action density to x^ is thermalized. This is clearly observed in SU{2), where a plateau 
at a value of x"^ exceeding P(/3) (perturbatively calculated) of the amount of the expected 
mixing, is detected ||^. The display of this plateau is possible in SU{2) since the average 
number of instantons in the first stage of the heating is nearly zero - instantons exhibit 
a very severe form of critical slowing down. Otherwise, they would start to give a 
contribution to the topological non-perturbative part of x^ before the thermalization 
of the mixing with the action density is reached, so preventing the observation of the 
plateau. 

It is interesting to investigate if this favourable situation happens also in the case 
of the 0(3) a model; here, the small-size divergence in the instanton size distribution 
is likely to cause a radical change of the scenario. Since the mixing term is negligible, 
a plateau is expected at a value of x'" corresponding to -P(/3). Previous results [Q 
0] seem to be in agreement with this expectation; however, the statistical fiuctations 
of those determinations prevent to absolutely exclude the eventual onset of instanton 
contribution during the first phase of the heating. 

In order to unmask the topological structure of the configurations in the heated 



ensemble, we use the "cooling" procedure |18|. It consists in a local minimization of the 



action with the purpose of destroying the quantum fluctuations, trying to preserve the 
background topological structure. We use an unconstrained cooling, consisting in the 
following replacement: 



0(x) -^ <p'{^) = c^y^ 
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-0(x + /i) - — 0(x + 2/i) 



(15) 



where a ensures the normalization of the new spin; this replacement exactly minimizes 
the action when the other spins are kept fixed. We measure the topological charge after 
each step of cooling. Typically, two behaviors are observed: the charge can rapidly go to 
zero, so revealing the absence of a topological background; otherwise, it reaches, after a 
short time, a constant value next to an integer, so indicating the presence of background 
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instantons. 

The irrelevant terms in the lattice action adopted - the Symanzik tree-level im- 
proved action - make the lattice action calculated on an instanton decrease when its 
size decreases. As a consequence, an instanton experiences under cooling a progressive 
shrinking, up to its destruction. The conclusion is that the cooling procedure affects in 
a certain degree the background topology, cutting off instantons with size smaller than 
a certain size, depending on the number of cooling steps performed. This number must 
be big enough to clear off the quantum noise, but not so much to completely erase most 
of the topological configurations. 

B. Heating an instanton 

Here, we apply the method suggested by Teper ||^ and already realized in [Q [Q , for 
the determination of Z{(3). Our contribution is an improvement of the statistic with 
the purpose to reveal eventual discrepancies from the two-loops perturbative calcula- 
tion; moreover, we study of the onset of small instantons over the starting topological 
background. 

We put an instanton of charge Qq = 1 and size 20 lattice units, in the middle 
of a 120 X 120 lattice. Starting from this configuration, we construct the ensembles 
{Ct} performing t sweeps of a standard heat-bath algorithm, and measure the average 
topological charge over each ensemble {Ct}. 

In Fig. 2 we present the results for three different values of j3. We observe that a 
plateau is reached after about 11 heating steps. The independence from (3 of the starting 
point of the plateaux has been observed also in previous works: it reveals that Z{I3) 
takes its origin from fluctuations of small size, whose thermalization does not undergo 
critical slowing down. 

In Table II we compare Z{j3) calculated at two loops with the non-perturbative 
estimate obtained by fitting values at the plateaux. We attribute the small discrepancy 
to further terms of the perturbative expansion of Z{[3): a fit gives zz = 0.097(8) and 



24 = -0.422(12), xV d.o.f. = 0.15. 

In order to check the extent of small instantons production during heating, we have 
analyzed a sample of 1000 configurations obtained after 15 heating steps at /5 = 1.45, 
where the small instantons contamination is expected to be the largest among the values 
of P considered. We classify each configuration assigning it to a definite topological 
sector according to its cooled charge. We find that ~ 30 configurations have left the 
Q = 1 sector: ~ 20 migrating into the zero charge sector, and ~ 10 into the Q = 2 
sector. This effect is to be attributed to the generation of small-size instantons (of 
charge —1 and +1 respectively), laying on the background topological configuration. 
We explain the small asymmetry observed as the effect of the explicit breaking of the 
charge symmetry by the starting configuration: starting from the Q = 1 sector, it is 
energetically more convenient to fall into the Q = sector than to rise into the Q = 2 
sector. We have verified that the systematic error induced on the determination of Z{l3) 
by a topological contamination of the observed extent is negligible within the statistical 
uncertainity of our result. 

C. Heating the flat configuration 

Here, in addition to the standard heat-bath algorithm, we use a faster local algorithm 
of thermalization, the overheat-bath 0. 

First of all, we have analyzed the thermalization process during heating. We have 
observed the behavior under heating of .^_l, the correlation length of the components of 
the spins orthogonal to the direction of polarization of the starting fiat configuration; 
with this choice we get rid of the effects of 0(3) symmetry breaking of the starting 
configuration. 



We have used the following definition |]T9 

1 



ei 



G±{0,0) 
4sin'(7r/L) [Gx(0,1) 

where 



(16) 



G±{k) = — ^(0x(a;) •0±(l/))exp 



}'i(.-y)-t 



(17) 



We have verified that this quantity reproduces at thermal equihbrium the value ^eq 
obtained with the standard definition. In Figs. 3, 4 we present some values of C,± 
versus t, the number of heating steps performed, obtained with the two algorithms 
of thermalization. We observe that the expected law C,±{t) = ct fits data nicely, 
with Ch = 0.335 ± 0.042, xV^-O-f. ~ 0.01, for the heat-bath, and Cov = 1-89 ± 0.16, 
X^/d.o.f. ~ 0.79, for the overheat-bath. We note that the overheat-bath is about six 
times faster than the heat-bath. 

In SU{2) Yang-Mills theory, topological fiuctuations are decoupled from the ordinary 
fluctuations as for thermalization time. This property allows the observation, during the 
heating of x^ starting from the fiat configuration, of an extended plateau corresponding 
to the non-topological contribution to x^y i-^- the mixing with the identity operator 
plus the mixing with the gluon condensate 0. Should this hold also in the 0(3) a 
model, a plateau in correspondence to -P(/3) would be detectable, and, as a result, a 
non-perturbative evaluation of such quantity could be obtained - the mixing term is 
here negligible. 

We have performed, at various values of (3, measures of x^ during the heating of 
the fiat configuration with the slower algorithm of thermalization, the standard heat- 
bath, each ensemble {Ct} containing from 10000 to 40000 configurations. Unfortunately, 
we could not clearly single out any plateau, always obtaining a steady drift to the 
equilibrium value (see Fig. 5). The situation is the same when the overheat-bath is 
used, with a six scale factor on the times of thermalization. 

Being aware of the ultraviolet nature of the model in study, the explanation of this 
behavior is straightforward: small-size instantons are precociously generated by the 
heating algorithm. 

In order to check this diagnosis, we sistematically cool ensembles {Ct}, obtained with 
a standard heat-bath, at values of t beyond the thermalization time of the perturbative 
fiuctuations indicated by the previous works on this matter [Q 0. We find that five 
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steps of uncostrained cooling are sufficient to unmask the topological content of the 
configurations through the observation of the behavior of the topological charge. If 
cooling is further protracted a part of instantons - the smaller ones - are destroyed: this 
phenomenon is revealed by the behavior of the charge under cooling that, after a short 
plateau at a value near to an integer, goes to zero in few steps. The result is that, as 
expected, the presence of topological fluctuations is relevant already at the ffist steps 
of the heating process: the argument that instantons are subject to a special kind of 
critical slowing down, more severe than non-topological fluctuations, seems not to apply 
to this model. 

Now, if topological configurations are sistematically removed from the statistical 
ensembles {Ct}, and x^ is averaged only over the trivial configurations, a long plateau 
should be observed after the thermalization of the perturbative modes. This plateau 



would correspond just to P(/3) as defined in Eq. (14). 

In order to get rid of instantons, we sistematically discard all configurations having, 
after 5 steps of unconstrained cooling, a topological charge greater than a fixed threshold 
value. An analysis at /? = 1.45 on a sample {C70} consisting of 100 configurations shows 
that this procedure eliminates almost all configurations with non-trivial topology. In 
Fig. 6 we show the results obtained at /5 = 1.45 for a chosen threshold: a clear long 
plateau is now observed. We have checked that data are stable under a small change of 
the threshold: this ensures that the plateau observed is a real effect and not an artifact 
of the subtraction procedure. We have repeated the above procedure for other values of 
P : a plateau is observed for t > t ^ 35, not depending on /? as expected; i is the time 
of thermalization of fluctuations with size / ~ 4 lattice spacings which, within our error 
bars, saturate the perturbative signal, P{(3)- We observe that P{(3), in comparison with 
Z{j3), is sensitive to larger size fluctuations: an indication of this could already be taken 
out from the perturbative series of -P(/3), which starts to get contribution from higher 
orders in 1//5. In Table III values of P{l3) fitted on the plateaux, Pnp{(3), are compared 
with those of PpertiP)- The latter quantity is obtained in Ref. [^] by a perturbative 
calculation combined with a fit on x^iP) ^it large (3: in this region, the probability 
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of occurrence of instantons is exponentially suppressed, as well the non-perturbative 
contribution of other origins, so PpertiP) is in all equivalent to a truncated perturbative 
determination of P{I3). The agreement observed in Table III reveals that P{I3) is a 
substantially perturbative quantity, well approximated by the first few perturbative 
terms. This result is not trivial, since in 0{N) a models perturbative series are not 



Borel-summable pO| 



Correlations between measures have been taken into account performing standard 
binning procedures; we observe in this regard that instantons, which are the main source 
of correlations for x^, have been rejected from the ensemble, and so the problem of 
correlations is less critical than usual. 

In order to study the approach to thermal equilibrium of x^ during the heating 
of the flat configuration, we exploit the faster overheat-bath algorithm. In Fig. 7 we 
show the behavior of x^ during the heating at /5 = 1.45: the time required to reach 
the equilibrium value is consistent with teg = ^Iq/cov This observation is confirmed, 
within our error bars, for other values of (3. The suggestion we get from this evidence 
is that instantons undergo the ordinary slowing down: the time of thermalization of all 
fluctuations depends only on their size - no matter of their topology. 

We further observe that the time x^ takes to overcome its perturbative value (see 
Fig. 8 and Table III) is the same for different values of (3. By a measure at /5 = 2.50, we 
have verified that this is just the time the algorithm takes to thermalize the perturbative 
fluctuations. The observed time (t ^ 8 heating steps) corresponds, according to the 
random walk law, to thermalized lengths of ~ 4 lattice spacings, in agreement with the 
results obtained with the heat-bath algorithm. The suggestion is that there is no non- 
perturbative contribution to x^ before perturbative modes are thermalized; otherwise, 
such non-perturbative contribution would give an additional boost to the signal, making 
it to overcome the perturbative threshold in advance at small /3 in comparison with the 
time taken at /9 = 2.50 to saturate the perturbative signal. 
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D. Small instantons size distribution 

On the basis of the results of the previous subsection, we make the following work 
hypothesis about the process of thermalization of topological fluctuations: we assume 
that the heating algorithm reproduces on the lattice at the time t the equilibrium size 
distribution of the average instanton number dN/dp, up to the length thermalized at 
that t, l(t)] we moreover assume that the quantity ^±(t), defined in (p^ , is a measure 
of l{t). In this hypothesis, the thermalized length obeys the law l{t) = \fci^ where c is 
the same constant that enters the random walk law for ^_|_ and that has been estimated 
for the two algorithms of thermalization in the previous subsection. Now, called iV(t) 
the average number of instantons in the ensemble at the time t, we write 

m)-l iP^- (18) 

In the hypothesis of non-interacting instantons, the average number of instantons 
in the statistical ensemble, A^, is simply related to the topological susceptibility: x = 
N/V. Indeed, denoting with n{n) the number of (anti-) instantons, and with Pn^ the 
probability to have a configuration with n instantons and n anti-instantons, the following 
relations hold 

X = ^J2(^-nfPn,n = ^J2(^ + n)Pn,n ^ ^ • (19) 

n,n n,n 

The hypothesis of non-interacting instantons plays a role in the second equality of (jT^) 

since, in this case, 

_ 1 {n + n)\ 
"•" - 2"+" n! n! "+" ^ ' 

(P/v is the probability of a configuration with A^ topological objects, instantons or anti- 
instantons). We note, by the way, that the above condition is satisfied if Pn,n = Pn Pn 
with Pn{n) poissonian, as in the semiclassical approximation. 

Eq. ([T9|) , when translated on the lattice in the framework of the field theoretical 
method, allows to obtain a measure of the average number of instantons in the ensemble 
{Ct}i N{t), from the measure of x^'- 
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N{t) xHP)-PiP) 



Anp 



So, we are in a position to check our hypotesis (|18|): in Fig. 9 we show N-^{t) as obtained 
from (^) at /5 = 1.45 together with the curve obtained from (^) with the choice 
dN/dp = A/p^'^^: we observe that the behavior of N{t) is nicely reproduced. 

The average number of instantons in the ensemble at time t can be directly obtained 
- without passing through x^ " from the analysis of the frequency of occurrence of 
configurations in the various topological sectors. Indeed, from the hypothesis of non- 
interacting instantons, it follows (see ([T9|)): 



N = T.Q'Pq ■ (22) 

Q 

In order to assign each configuration to its topological sector Q, we exploit the above 
described cooling procedure. 

In Fig. 10 we show the behavior of Nqlt) so obtained at /? = 1.45. We see that NQ{t) 
is ~ for t < 2 heating steps, corresponding to instantons of size ~ 2 lattice spacings: 
instantons smaller than this size do not survive the cooling process. 
For 2 < t < t (t is the time of thermalization of perturbative fluctuations, corre- 
sponding to ~ 4 lattice spacings), NQ{t) is different from zero even if x^ has yet no 
non-perturbative contribution, as observed in the previous subsection. We conclude 
that instantons with size lower than ~ 4 lattice spacings - though present in the ensem- 
bles {Ct} - give no non-perturbative contribution to x^^ so having the same effect of 
perturbative fluctuations. 

For t > i, NQ{t) is systematically greater than N^{t): this discrepancy is simply ex- 
plained, since x^ loses the part of the topological contribution coming from instantons 
with sizes ranging from ~ 2 to ~ 4 lattice spacings. In formulae: 

/■'(*) dN 



while 



NQ{t) = 1^^ dp ^ . (24) 
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The mismatch is expected to be just 

■^4 dN 

dp— ^ NQ{i) , (25) 

2 dp 

for all t > i. We find that the quantity NQ{t) — Ny^{t) is well fitted by a constant 
(x^/d.o.f ~ 0.1) whose value is consistent with the prevision ( p5D - see Fig. 11. The 
same happens at /5 = 1.525. 

On the basis of the assumption (^) on the mechanism of thermalization of instantons 
during heating, we can extract from the lattice the size distribution of instantons in the 
physical vacuum, dN/dp. Indeed, from Eq. ([ISD it follows 

dN 
dp 



1 /dN\ 
^ = ,(,) " dl{t)/dt [H) 



In Fig. 12 we plot dNq/dp and dN^^/dp versus p aX (3 = 1.45 as obtained from Eq. (^61) . 
We find that the two determinations of dN/dp nicely overlap. The same is done in 
Fig. 13 at /5 = 1.525; data at /? = 1.525 have a suppression oc a^ in comparison with 
data at /3 = 1.45, since dN/Vdp is a renormalization group invariant quantity when V 
and p are measured in physical units. 

In Fig. 14 we collect all our data for dN/Vdp expressing all quantities in physical 
units. Data coming from different values of (3, covering a spectrum of p ranging from 
0.2 to 0.9 units of the inverse mass gap, dispose within statistical fiuctuations along the 
same curve, as expected. A global fit with a function A/p"-' shows that n from 1 to 2 is 
consistent with data. In fact, n seems to increase with p: initial data give n ~ 1.5, while 
data at p > 0.4 give n ~ 2. This behavior is intermediate to the expectation n = 1, from 
semiclassical prediction at small sizes, and the result n = 3 recently obtained in [pT[], 
where dN/Vdp is determined with different techniques in the region 1 < p < 5.5. These 
results seem to suggest for dN/Vdp an ever more severe suppression as p increases. 

Up to the sizes investigated by our analysis (p ~ 0.2 inverse mass gap units), no 
physical ultraviolet cut-off has been observed in the instanton size distribution. Due to 
this ultraviolet dominance, Xnp loses a noticeable part of the non-perturbative contri- 
bution, if the discretization is too coarse: we argue that this is in fact the case for the 
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values of (3 of the numerical simulations up to now performed. Since the amount of the 
lost contribution decreases as 13 increases, a scaling defect for Xnp should be observed. In 
the following, we show the results of Xnp/ PiP) ^^ thermal equilibrium for three values 
of (3 in the expected scaling region; /(/3) = 2'Kj3e~'^'^^ is the two-loop renormalization 
group function and Xnp ^^^ been obtained from standard Monte Carlo simulations using 
the non-perturbative determinations of P(/3) and Z{f3): 

(5 = 1.45 : 196 ±11 (27) 

13= 1.50 : 224 ±13 
p= 1.525 : 250 ±14 . 

The increasing behavior observed indicates that, for the values of j3 in study, Xnp 
has not yet saturated all topological contributions. 



V. CONCLUSIONS 

In this paper we have studied the behavior of topological small-size fluctuations un- 
der heating in the 0(3) a model. The main conclusion is that small instantons undergo 
ordinary slowing down: instantons with a given size reach the thermodinamical equi- 
librium distribution when the heating algorithm attains the corresponding length. A 
different conclusion in this regard is settled in SU{2) Yang-Mills theory, where instan- 
tons are believed to undergo a form of critical slowing down more severe than ordinary 
fluctuations. We think that this state of things is to be attributed to the pathological 
characteristic of the topology in the 0(3) a model, exhibiting ultraviolet dominance. 

As a by-product of the present work, we have obtained a high-precision non- 
perturbative estimate of the mixing of the lattice topological susceptibility with the 
unity operator, flnding good agreement with perturbative results. This shows that this 
quantity is essentially perturbative. We have also obtained the size distribution function 
of the average number of instantons in the physical vacuum in a region of sizes not yet 
explored, reaching sizes of a small fraction of the inverse mass gap, the physical length 

15 



of the model. These results lay on a two-fold assumption: first, instantons thermalize as 
ordinary fluctuations; second, instantons do not interact among themselves. The first 
assumption is consistent with other results of the present work; the second is reasonable 
if one takes into account that the typical sizes of instantons involved in our analysis are 
much smaller than the lattice size. 

Up to the sizes explored, no physical cut in the distribution function has been ob- 
served, so having a new evidence that the topology of the 0(3) a model is strongly 
ultraviolet dominated. Moreover, we have found that the lattice-regularized version of 
the topological susceptibility has a non-physical cut-off in the size of instantons cor- 
responding to ~ 4 lattice spacings. As a consequence, a part of the non-perturbative 
contribution to the topological susceptibility is lost on the lattice, at least in the region 
of 13 up to now explored in numerical simulations. We attribute to this loss of non- 
perturbative signal - decreasing as (3 increases - the observed deviation from scaling of 
the topological susceptibility on the lattice. 

We have also argued on the basis of a perturbative calculation, that the mixing of 
the lattice topological susceptibility with the action density is a negligible part of the 
total non-perturbative signal. 

The ultimate consideration of this work is a caveat: when the 0(3) a model is used 
as a laboratory to get insight into methods conceived for the theory of physical interest, 
i.e. QCD, the possibility to be wrecked in the unphysical artefacts of the model must 
be seriously considered. 
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APPENDIX 
Calculation of the mixing coefficient 

For both continuum and lattice calculations at low temperatures, the perturbative 



expansion is obtained by setting (p = [ttj, Wl — Djvr^ ], i = 1, 2. The perturbation 
theory suffers from infrared divergences, which can be cured by adding a magnetic term 
to the action 



Sm = Jd^xh ll-Y^nf . (28) 

Indeed, Sm explicitly breaks the 0(3) invariance and acts as a mass term for the vr 
field. The 0(3) invariant quantities (and the relations among them) are free of infrared 
divergences, and have a well-defined limit for /i — > 0. 

From the standard theory of the renormalization of composite operators it follows 
that the coefficients Z{f3) and A{f3) enter the relation 

rf^ nLf ,nL, ,1 = a'zipf rj';)^'^^) ,, + aM(/3) rS ., ,^ + 0(0^) , (29) 

[2_^^Q^{x)Q^{y%Ten V/^^ [j^Q{x)Q{y)],ren ^^^ [S(y)\,ren \ I ■: \ J 

where r^^j ^^^ stands for n-point {n ^ 0) renormalized F-function with the insertion 



of the operator O and MS is an arbitrarily chosen renormalization scheme. Eq. ([29| ) 
allows to calculate A{(3), order by order in perturbation theory, as a series in 1/(3: 
-A{(3) = Sn'^n//^"- I^ Eq. (|25|), an integration over the variable y is allowed since there 
are no operators other than S{x) with the same quantum numbers of x (this would be 
forbidden in the case of mixings with several operators differing by an integration by 
parts). After the integration, Eq. (ESI) becomes 



^mWUen = ^(/3)^ r(5)J^^),„ + A(/3)F[^)„),__ + 0(a^) , (30) 

where V = L^a"^ is the space-time volume. Now the vertices of the topological charge 
appear in the diagrams at zero momentum, with a considerable simplification of the 
calculation. 
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In view of a perturbative evaluation of A{f3), the most suitable choice is n = 2, with 
an insertion of two external lines of the field vr with the same index i. At the lowest 
order: 

r[?(0)],.en = V , (31) 

where p is the momentum flowing in the external lines, and there is only one diagram V - 
order 1/(3^ - contributing to rLa /yi, both on the lattice and on the continuum (Fig. 1). A 
straightforward calculation shows that the diagram vanishes on the continuum: V\cont = 
0. As a consequence: 



Vliatt = -Jr2p' + 0{a')- (32) 



a^ is just the coefficient we are looking for. 

We now present the complete calculation of V\iatt- All the results will be expressed 
as obtained in the infinite-size lattice: they can be translated in the finite-size case by 
simply interpreting integrals as discrete sums. We first give the explicit expression of 
the four-legs vertex of the lattice topological charge operator Q^: 

^/n.^ V{p^ , p2) {2nf5i ^p, ) 7ri(pi)7r2(p2) 7r{ps) ■ n{p,) , (33) 

where all the integrations are from — tt to vr, the index i runs from 1 to 4, and 

y{Pi^P2) = f^fMu [smpif,smp2^ - sin{pi+p2)^{smpi^ -sinp2u)]. (34) 

The result of the diagram in Fig. 1 is: 

1 f (fpi (Pp2 lip , Pi ,P2) 



(35) 



87r2/53 J {2ny (27r)2 (D,^ + h){n^, + h){n^+,,+^, + h) ' 
where Dp is a shorthand for the inverse propagator of the Symanzyk tree-level improved 



action 



°p = E^i + -^UPir , (36) 



Pi = Pf.Pl , Pm = (^'"'^ - 1 
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and I{j) ,pi ,p2) is a compact notation for 

V{p,pi)V{p,pi) + V{p,pi)V{p,p2) + 2V{p,p2)V{pi,p2) + 

{V{pi,p2)f + V{pi,p2)V{p + pi+p2,-p2)- (37) 

It may be checked that the integral in (^) is convergent in the hmit /i — > 0. The 
correct procedure would be now calculating exactly the above integral and extracting 
the term proportional to p^ in the result of the integration. A look at the complete 
expression suffices to discourage any attempt in this direction. The other, apparently 
viable, possibility of performing Taylor expansions in the integrand and isolating the 
p"^ term before integration, does not work [Q. We have found a way to get the result 
without any Taylor expansion, relying only upon the symmetries of the vertex V{pi,p2). 
Let us consider, for example, the contribution to ( p5|) coming from the third term in 



m 



V{p,P2)V{pi,p2) _ 1 ^ ,gg. 



After the translation pi ^ pi — p we have 

X = 2 f ^(P^ P2)V{P1-P, P2) .gg^ 

We observe that V^(p, P2) factors out a sinp ~ p term, while V{pi — p , P2) may be 
decomposed into V{j)i,p2) and the remaining part, ^-(^,^1,^2)7 which is 0{j)). When 
^{Pi'P2) multiplies V^ the p^ term comes out trivially and it is sufficient to put p = 
in the propagators; the part with ^^(^,^2)^(^1,^2) reduces to — X3 after the following 
transformation of the integration variable pi. pi —>■ pi — P2, Pi —>■ — Pi and using the 
property V{-pi - P2 , P2) = - V{pi,p2). 

A similar procedure works for the contribution to (|35|) from the fourth and the fifth 
terms in (^), while the first and the second terms need no manipulation because the 
factorization of p"^ is trivial. 

The sum of the various terms in (^) after the factorization of 2p^ is an IR-convergent 
integral where h can be safely put equal to zero; the mixing coefficient a^ we are looking 
for, is: 
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it^I °Pi°P2°Pi+P2 {sin^PiM(2 + cospi^)2 + sinpi^sinp2/.(2 + cospi^)(2 +cosp2i.) 
sinp2i.(2 + cosp2m)[ (sinpi^, - smp2u) cos(pi +P2)^^ - {smp2u - sin(pi + p2 )u) cospi^ ] 
+i [ (sinpi^ - smp2u) cos(pi +^2)^ - (smp2i/ - sm{pi +P2),,) cospi^ J^j , (40) 

where there is no sum over repeated indices; /i and z/ indicate two different directions. 

In Table I we report values of 03 for various lattices. 

In order to evaluate the relevance of the mixing of x^ with the action density in the 
determination of x from lattice simulations, we now need an estimate of the condensate 



(5(0)). The only available result comes from the large A^ limit IT?]: 

{SmiN^oo = -^m' . (41) 

By using the exact relation between the mass gap of the theory, m, and the scale 
invariant parameter of the lattice Symanzik- improved theory, Asy [H], we obtain the 
estimate for A^ = 3 

(S(0)) ~ -200 A|y . (42) 

This result, combined with our estimate for A{f3), gives a mixing term which is negligible 
with respect to the topological non-perturbative signal: a.t (3 = 1.45, for instance, 
A{(3){S{0))/[Z^{f3)x] ~ 6 • 10"^. We believe that higher order corrections in 1/A^ and 
1/(3 cannot spoil the validity of this argument. 
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TABLE I: 03 for a variety of lattices. 03 is the first non-vanishing term in the perturba- 
tive expansion A{P) = J2n (^n/P^ of the mixing coefficient of the topological susceptibility 
with the action density. L is the lattice size. 



L 


as X 10^ 


L 


as X 10^ 


L 


as X 10^ 


10 


0.31450626 


60 


0.33496730 


110 


0.33538507 


20 


0.33022955 


70 


0.33512505 


120 


0.33541334 


30 


0.33318636 


80 


0.33522745 


130 


0.33543534 


40 


0.33422462 


90 


0.33529767 


140 


0.33545280 


50 


0.33470577 


100 


0.33534790 


00 


0.33556212 



Table II: Z{j3) versus [3. Z{(3)2ioops is the multiplicative renormalization calculated to 
two loops. Z{j3)mc is the multiplicative renormalization calculated by heating an in- 
stanton; the size of the lattice is 120; Stat is the statistic of the simulation. Since data 
on the plateau are correlated, we report the average error of data. 



p 


Z{(3)2loops 


Z{P)mc 


Stat 


2. 


0.643 


0.6285(26) 


16000 


1.65 


0.564 


0.5298(42) 


16000 


1.45 


0.500 


0.4362(31) 


64000 
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Table III: P{f3) versus [3 for different choices of the threshold. P{f3)pert includes the 
perturbative tail calculated to four loops and higher order terms fitted from thermal 
equilibrium Monte Carlo data at large j3. P{j3)np is the perturbative tail calculated by 
Monte Carlo techniques; the size of the lattice is 120; Stat is the statistic of the simu- 
lation. 



/? 


Threshold 


P{P)pert X 10^ 


P{P)np X 10^ 


Stat 




0.2 




5.24(6) 




1.45 


0.4 
0.6 


5.21(14) 


5.26(5) 
5.28(6) 


5000 




0.2 




4.22(3) 




1.50 


0.4 
0.6 


4.20(12) 


4.23(3) 
4.24(3) 


5000 




0.2 




3.82(2) 




1.525 


0.4 
0.6 


3.89(10) 


3.83(2) 
3.86(2) 


5000 




0.2 




3.50(2) 




1.55 


0.4 
0.6 


3.54(9) 


3.51(2) 
3.53(2) 


5000 
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FIGURES 



FIG. 1. Diagram contributing to the two-point proper function with insertion of 
{Q^Y /V . White crossed blobs indicate the charge operator. 

FIG. 2. Q^ versus the heating step. The thermalization is performed by the heat- 
bath algorithm aX (3 = 1.45, 1.65, 2.00 on a 120^ lattice. The solid lines indicate the 
values of Z{f3) estimated by averaging data at the plateaux. 

FIG. 3. ^\ versus the heating step with the heat-bath algorithm dX j3 = 1.45. The 
solid line indicates the best fit with the law ^^ = ct (with the relative error). 

FIG. 4. As in Fig. 5 with the overheat-bath algorithm. Data are obtained at two 
different values oi 13: 1.45 and 1.525. 

FIG. 5. x^ versus the heating step at /3 = 1.45 starting from the flat conflguration 
on a 120^ lattice. The thermalization is performed by the heat-bath algorithm. The 
two solid lines indicate the estimate of P{I3) by perturbative calculation and the value 
of x^ obtained at the thermal equilibium (with the respective errors). 

FIG. 6. x^ versus the heating step at /9 = 1.45 starting from the flat conflguration 
on a 120^ lattice. The thermalization is performed by the heat-bath algorithm with 
the procedure of subtraction of topological conflgurations described in the text. The 
threshold is 0.4. The solid line indicates the value of P{l3) estimated averaging data at 
the plateau (with the relative error). 

FIG. 7. The approach to equilibrium of x"^ at /? = 1.45 with the overheat-bath 
algorithm. The two solid lines are the non-perturbative estimate of P{I3) and the value of 
X^ obtained by standard MC simulations at the thermal equilibrium (with the respective 
errors) . 

FIG. 8. The overcoming of the perturbative value during the heating of x^ at various 
j3. The solid lines are the perturbative values. 

FIG. 9. The average number of instantons obtained from x^ at /3 = 1.45 versus 
the heating step. The solid line represents the predicted behavior in our model of 
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thermalization (see the text). 

FIG. 10. The same quantity of Fig. 9 obtained from the frequency of the configura- 
tions in the various topological sectors. 

FIG. 11. Comparison of the quantities of Figs. 9 and 10 after the subtraction from 
data of Fig. 10 of the instantons with sizes smaller than the the typical perturbative 
ranges, which give no contribution to the non-perturbative part of x^ ■ 

FIG. 12. The instanton size distribution (in lattice units) as obtained from the two 
methods described in the text, aX jS = 1.45. 

FIG. 13. As in Fig. 12 with ^ = 1.525. 

FIG. 14. The instanton size distribution in physical units obtained collecting data 
at different values of (3. The solid line is the best fit with n = 1.5. 
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